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Abstract 

The Josephson current through a ID quantum wire with Rashba spin-orbit and electron-electron 
interactions is calculated. We show that the interplay of Rashba and Zeeman interactions gives 
rise to a supercurrent through the ID conductor that is anomalous in the sense that it persists in 
the absence of any phase difference between the two superconducting leads to which it is attached. 
The electron dispersion asymmetry induced by the Rashba interaction in a Luttinger-liquid wire 
plays a significant role for poorly transmitting junctions. It is shown that for a weak or moderate 
electron-electron interaction the spectrum of plasmonic modes confined to the normal part of the 
junction becomes quasi-random in the presence of dispersion asymmetry. 
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I. INTRODUCTION 



In recent years the concept of a Luttinger liquid (LL) as a realistic model of interacting 
electrons in one-dimensional (ID) metallic structures has received experimental support 
(see e.g. Refs. . Quantum wires (QW) in laterally constrained 2D electron gases 

(2DEG) and single wall carbon nanotubes (SWNT) - are the two best known structures 
where LL behavior has been established both theoretically and experimentally. 

In SWNTs, where the interaction effects have been shown to be stron g 3 i 4 i 5 , interactions in- 
fluence the charge and spin transport through the nanotube. When a repulsively interacting 
LL is coupled to leads (M) of noninteracting electrons, two qualitatively different regimes of 
charge transport may be realized depending on the quality of the LL/M electrical contacts. 
For tunnel contacts charge transport through the system is strongly suppressed at low tem- 
peratures and bias voltages^ by the repulsive electron-electron (e-e) interaction. In contrast, 
for adiabatic contacts when electron backscattering is negligibly small, the conductance is 
not renormalized by the interaction^^. 

These two types of charge transport behavior also characterize the superconducting prop- 
erties of a LL wire coupled to superconductors. For adiabatic contacts only Andreev scat- 
tering of electrons occurs at the boundaries between LL and bulk superconductors (LL/S 
boundaries). This process does not lead to a redistribution of charge density along the wire 
and therefore the Coulomb interaction does not influence the supercurrent through a perfect 
LL. The above result was proved in Ref. Q| by a direct calculation of the Josephson current 
through a long S/LL/S junction, L ^> £ , (L is the junction length, £ — %Vf/\&\ is the 
superconducting coherence length, A is the superconducting order parameter). In the case 
of tunnel S/I/LL/I/S junction — where "I" denotes the insulating "layer" — the repulsive 
e-e interaction results in a renormalization of the junction transparency and the critical 
Josephson current is strongly suppressed^. 

Here we consider the influence of spin-orbit (s-o) interaction on the Josephson current 
through a long S/I/LL/I/S junction. It has been known for a long time that the s-o inter- 
action is strong in a 2DEG formed in a GaAS/AlGaAs inversion layer (the Rashba effect) 
and that it can be controlled by a gate voltag o 13 i 14 i 15 . In what follows we will consider a 
quantum wire in a laterally confined 2DEG coupled to superconducting electrodes via tunnel 
barriers. 
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FIG. 1: Schematic energy spectrum of ID spin-1/2 electrons with dispersion asymmetry. The 
subbands "1" and "2" are characterized by their Fermi velocities Vip ^ v 2 f- in the case of weak 
Rashba interaction spin projections in the subbands for each given momentum are opposite. For 
strong Rashba interaction spins in subbands are parallel and they are opposite for right- and 
left-moving particles. 



The influence of the Rashba effect on the electron spectrum and on the 



transport prop- 



16 



17| . where it 



erties of quasi- ID quantum wires has been studied theoretically in Refs. 
was shown that the s-o interaction not only splits the electron spectrum into spin-" up" 
and spin-" down" subbands, but additionally breaks the chiral symmetry. This implies that 
left- and right-moving electrons with the same spin projection have different Fermi veloc- 
ities. Since the time invariance (T-symmetry) of the spin-orbit Hamiltonian implies that 

( ( F) ( F*\ ( F*\ 

v r\ = v li = v if and v R!i = v L ^' = v 2 p the strength of the Rashba effect in a single- channel 
QW can be characterized by a dispersion asymmetry parameter A a = {yip — v 2F ) / {vip +v 2F ) . 
In Refs. jlfi|l7j | it was assumed that in the presence of Rashba interactions the electron spins 
in a quasi-lD wire are aligned as in the 2D case (see FigGJ solid lines for spin projections). 
Although this assumption is not valid for a strong Rashba coupling^, the model considered 
by the authors of these references is interesting in itself. It allows one to study the effects 
of dispersion asymmetry (A a 7^ 0) on the electron dynamics and in the limit A a — ^ it 
reproduces the standard results for spin-1/2 electrons without s-o interaction. 

Since the electron spin is not conserved in the presence of s-o interactions the classification 
of spin states assumed in Refs. jlfil|l7j | is not evidently correct. Actually, as was shown in 
Ref. 18] , it can be justified only for a weak Rashba interaction. In the most interesting case of 
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strong Rashba interaction, when the characteristic energy scale introduced by s-o coupling is 
compared with the energy spacing of ID subbands, the average spin projections for electrons 
with large (Fermi) momentum are different. The total energy is minimized when all right- 
moving (R) electrons have parallel spins and they are in the opposite direction to the spins 
of left-moving (L) electrons^. In what follows we choose the sign of the Rasha interaction so 
that R-electrons have "down"-spin and L-electrons are "up"-spin particles (see FigHJ dashed 
lines for spin projections). Notice that under conditions when the Rashba effect is active 
the electron spin lies in a (2D) plane and orthogonal to the electron momentum (in the ID 
case this direction is fixed and "up" and "down" spin projections are well defined). 

At first we consider the influence of electron dispersion asymmetry in the model elabo- 



rated in Refs. 
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171 ] on the superconducting properties of S/I/LL/I/S junction. We calculate 



Josephson current perturbatively on the junction transparency D = \tit r \ 2 , (|t«, r | 2 <C 1 are 
the transparencies of tunnel barriers at left and right LL/S interfaces) and for arbitrary 
strength of electron-electron interaction, dispersion asymmetry A a and Zeeman splitting 
Az = g^sB (g is g-factor, /xb is Bohr magneton and B is the magnetic field). Two different 
geometries of S/LL/S junction are considered. In the first case an effectively infinite LL is 
connected by the side electrodes to the bulk superconductors (Fig|2J). In this geometry^ 
one can use periodic boundary conditions for plasmons and all calculations can be done 
analytically even in the presence of s-o interaction. When dispersion asymmetry is negli- 
gibly small (A a — > 0) we reproduce the formula for Josephson current derived in Ref. 11 1. 
As was shown in the cited paper the influence of Coulomb interaction on a supercurrent 
through a tunnel junction J = J c simp results in suppression of the critical current, which 
in the presence of Zeeman splitting takes the form J c = J^°>R i (g c ) cos(A^/A I/ ), where 
DAl (Al = Hvf/L) is the critical current for noninteracting electrons and the 
interaction- induced renormalization factor Ri (the subindex "i" labels the case of effectively 
infinite LL wire) is small for repulsively interacting electrons Ri(g c <C 1) <C 1 (here g c is 
the LL correlation parameter in the charge sector). The Zeeman interaction in the absence 
of Rashba effect results only in additional sign alternating harmonic factor in the critical 
current (see also Refs. 

We assume that the magnetic field is local and influencing only the normal (nonsuper- 
conducting) part of the junction. (This can be realized in an experiment for instance with 
the help of a magnetic tip in a scanning tunnelling microscope). The interplay of dispersion 
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FIG. 2: SLLS-junction of lengh L formed by an effectively infinite Luttinger liquid coupled to bulk 
superconductors by side electrodes. 

asymmetry and the Zeeman interaction leads to beatings in the supercurrent considered as 
a function of the local magnetic field B. A more unusual prediction is the appearance of 
supercurrent J a even at ip = 0. The existence of this anomalous Josephson current is related 
to the breaking of chiral invariance in quasi- ID quantum wires and the effect manifests itself 



already for noninteracting particles (see Ref. 2l|). 



A more realistic geometry for an S/LL/S junction is a finite LL wire (of the length L) 
coupled via tunnel barriers to bulk superconductors (FigJSJ). We will assume that the barrier 
transparencies are unequal and small (nonsymmetric tunnel junction \ti\ ^ \t r \ <C 1) and 
evaluate the (^-dependent part of the ground state energy in perturbation theory using the 
junction transparency D = \tit r \ 2 as expansion parameter. Normal and Andreev scattering 
at the interfaces can be taken into account by the boundary terms^ in the Hamiltonian of 
the S/I/LL/I/S junction. To first order in the junction transparency the problem is reduced 
to the evaluation of four-fermion correlation functions for a two channel LL Hamiltonian 
with the boundary conditions implying the absence of particle current through the S/LL 
interfaces at x = 0, L. In the absence of spin-orbit interaction the problem of quantization 
of plasmon modes in a finite LL with open ends was solved in Ref. |22j. Here, we generalize 
the quantization procedure proposed in the cited paper to the case of spin- 1/2 fermions with 
dispersion asymmetry. 

The two-channel LL Hamiltonian describing our system is diagonalized exactly by the 
canonical transformation suggested in Ref. We show that the spectrum of plasmons in a 
LL with open ends in the presence of dispersion asymmetry is determined by a transcendental 
equation. In the general case the spectrum forms a set of quasi-random energy levels. For 
noninteracting electrons, or when the energy dispersions are symmetric (vxf — V 2F — vf), 
the spectrum is reduced to a set of equidistant energy levels. In the limit of strongly 
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FIG. 3: A Luttinger liquid wire of length L coupled to bulk superconductors via tunnel barriers 
with transparencies Du r y 

interacting particles the plasmon spectrum also becomes regular. We calculate the Josephson 
current in the cases when the spectral equation can be solved analytically. 

We find that for noninteracting electrons the critical Josephson current through a tunnel 
S/QW/S junction is enhanced by the presence of dispersion asymmetry. This behavior is 
specific for ID electrons and the effect disappears in 2D junctions^. 

When the Rashba effect is not pronounced (A a — ► 0), the Josephson current through a 
tunnel S/LL/S junction is described by an expression analogous to the formula derived for 
an effectively infinite LL. However, the interaction- induced renormalization (suppression) of 
the critical supercurrent is much stronger for a finite LL than for an infinite one Rf{g c "C 
1) <C Ri(g c "C 1) provided the electron-electron interaction is short-ranged. 

As already mentioned in the Introduction, the electron spin states in a ID quantum 
wire in the regime of strong Rashba effect are fully determined by s-o interaction and the 
electrons with large (Fermi) momenta behave as truly chiral particles. That is the electron 
spin polarizations and the direction of their motion (right/left) are strongly correlated and 
all right (left)- moving particles, irrespective of their Fermi velocities, have parallel spins^ 
which are opposite to the spin polarizations of left(right)-moving electrons. So it is reason- 
able to expect that in this case the magnetic field via the Zeeman interaction will induce 
an anomalous supercurrent (at (p — 0) even in the absence of dispersion asymmetry. We 
calculate the Josephson current in a SILLIS junction in a model when the Rashba s-o inter- 
action is smoothly switched on in a ID QW and spin-flips are not accompanied by electron 
backscattering. An anomalous influence of Zeeman splitting on the critical supercurrent is 
predicted. 



6 



II. PROXIMITY-INDUCED SUPERCONDUCTIVITY IN A LUTTINGER LIQ- 
UID WIRE WITH CHIRAL SYMMETRY BREAKING 

It is physically evident that the Coulomb interaction in a long S/I/LL/I/S junction 
suppresses the critical supercurrent due to a strong Kane-Fisher renormalization of the bare 
tunneling matrix elements. The Josephson current through a Luttinger liquid coupled to 
bulk superconductors via tunnel contacts was first calculated by Fazio et alM who showed 
that the critical supercurrent is multiplicatively renormalized (suppressed) by a repulsive 
electron-electron interaction. The calculations were performed in linear (FigEJ and ring-like 
geometries. In both cases periodic boundary conditions for the plasmonic modes can be 
imposed. Although from an experimental point of view the considered geometries of an SNS 
junction look rather artificial, they do allow one to simplify the calculations. 

For noninteracting electrons the critical supercurrents in an SNS junction formed by a 
long (effectively infinite) quantum wire connected to superconductors by side tunnel contacts 
(separated by a distance L) and in an SNS junction where a finite length QW bridges the gap 
(of the same length L) between two superconductors differ only by a numerical factor. If the 
QW is treated as a Luttinger liquid this factor becomes a function of the interaction strength 
and can be evaluated analytically (see below). When both electron-electron interactions and 
dispersion asymmetry are present the calculations are more cumbersome. We start with the 
case of a side-contacted QW where we are able to analytically evaluate the supercurrent for 
arbitrary interaction strength and dispersion asymmetry parameter. 

The Hamiltonian H = Hll + Hf, of a S/I/LL/I/S junction is a sum of the LL Hamilto- 
nian Hll and the boundary Hamiltonian H^. The latter describes the effective boundary 
pairing and scattering interactions produced by the superconducting and normal scattering 
potentials at the points x = and x = L (see Ref. |26|). When chiral symmetry is broken 
the corresponding spin- 1/2 LL Hamiltonian expressed in terms of charge densities of chiral 
fields takes the form 



+ —c(PR\PRX + Pl]Pli + PriPl] + PriPli 

+ PR\PLl+ PRIPL\)} ■ 

where Pr/l,i/i are the charge density operators of right /left-moving electrons with up/down- 




(1) 
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spin projection, Vo is the strength of electron-electron interaction (Vo ~ e 2 ) and «i( 2 ) = 
v i(2)F + Vo/2nh. The Fermi velocities V\p ^ v 2 f are different due to an assumed electron 
dispersion asymmetry (see FigHJ. We have neglected the magnetic field- induced corrections 
to the Fermi velocities and assumed that the effective electron-electron interaction has no 
significant magnetic field dependence. Both the neglected effects are of "1/ep" -order (see 
e.g. Ref. Q) 

and they are irrelevant for Zeeman splittings A z < |A| < e F . 
The Hamiltonian (HI) is equivalent to a two-channel LL Hamiltonian and can be di- 
agonalized by the canonical transformation suggested in Ref. |23| (see Appendix I). The 
diagonalized Hamiltonian is 

H d = nhJ dx{si{p 2 Rl + p 2 L1 ) + s 2 {p 2 m + p 2 L2 )}, (2) 

where S(i,2) are the velocities of noninteracting bosonic modes (see Appendix I). 

We assume strong normal backscattering at the S/N-boundaries (tunnel junction). In 
this limit the pairing Hamiltonian contains a small factor - the amplitude of Andreev 
backscatterin g 28 1 29 

rj'° ~ D ril exp + (f r ,i) , (3) 



1 U + ' 

where -D r (z) C 1 is the transparency of the barrier at the right (left) interface, (p r m is the 
phase of the superconducting order parameter on the right (left) bank of the junction. The 
boundary Hamiltonian for our two-channel system can be expressed in terms of the Andreev 
scattering amplitudes Eq. © up to an overall numerical factor C, which will be specified 
later 

H b /C = hv 1F [rfy m (0)V Ll (0) - rf^ m (L)^ Ll (L)] (4) 
+ hv 2F [r* (0 * fli (0)* iT (0) -r* (r) * m (L)* LT (L)] +h.c. 

To second order in the Andreev scattering amplitude the phase dependent part of the ground 
state energy takes the form 

6E p ){(p) = g I^M' 2 = 1 jf° dr(0\Ht(r)H b (0)\0) , (5) 

where H b {r) is the boundary Hamiltonian (J3J in the imaginary time Heisenberg represen- 
tation. After substituting Eq.(J3J) into Eq.fjSJ) we get the following expression for 5E^ 2 '((p) 
expressed in terms of electron correlation functions 

( POO 

5E^ 2 \ V ) = -4ChM\rfr$> jf dr[v 2 F (^ m (r, Q)9 Ll (r, 0)^(0, L)¥ m (0, Q) + «T«H>] 

(6) 



We will calculate the electron correlation functions in Eq. © by making use of the 
bosonization technique. Notice that the Zeeman splitting introduces an extra x-dependent 
phase factor in the chiral components of the fermion fields. This interaction can be taken 
into account (see e.g. Ref. by replacing the fermion operators in Eq. (jUJ) by ^1^, where 

9$ = e*p{iK z z)9„ , K z = ^^ t (7) 

Here v F = (v 1F + v 2 f)/2, fi = (R, L) = (1, -1) , a = (|, 1) = (1,-1) , A z is the Zeeman 
splitting, and A a = {v\p — v 2 f) / {v\f + V2f) is the parameter which characterizes the strength 
of chiral symmetry breaking. 

The standard bosonization formulae now read 



xV f x ex P {±iV4n$ R{m (x,t)} exp{±iy/±*$ R(L)d (x,t)} 

^/27rai( 2 ) ■ s J2Tca 2 (i) 

(8) 

where a% }2 are the cutoff parameters of the two-channel LL. The chiral bosonic fields in 
Eq. (JHJ) for a finite length LL are represented as follows (see e.g. Ref. j^ll ]) 

$*(£),tOM) = ^fi(L),T + H T X T Vl{2)t + <Pr(l),i {x, t) , (9) 

2 M(2) 

*RtLu(x,t) = UniLU + V ^ + <PR(L U (x, t) ■ (10) 

2 -^2(1) 

Here the zero mode operators ((pR(L),a, ft^) obey the commutation relations [^mo-, fl a i] = 
~^i& a ,a' an d the nontopological (harmonic) components (fR(L),j(x,t) are 



<PR(L)j[x 



^)-T.^M q{±X - Vit \ + h.c) , (11) 

where & 9 (&J) are the standard bosonic annihilation(creation) operators. The effective quan- 
tization lengths Lj (j = 1, 2) depend on the boundary conditions and will be specified in the 
next section. 

As is well known (see e.g. Ref. j^), the topological excitations for an effectively infinite 
LL play no role and can be omitted in Eqs. © and (fTU|) . After strightforward transformations 
Eq. (jHj) is reduced to the following expression 



{A„ roc A „ roc 

v\ F cos(^ - - A / drill (r) + v\ F cos(<^ + — ^-) / tMI 2 (r) 
Ai L -/o " ' A 2 l ->o 



(12) 
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where D = D[D T is the junction transparency, Ai( 2 )£ = Uv^f/L and 

ni (2) (r) = 1 exp{27r[((y g (r,-L)y g )) + ((e p (r,-L)e p ))±((e p (r,-L)y g )) 
(27rai( 2 )) 2 

± «^(t, -L)6 p »]} . (13) 

Here <^ CT = ^(O, 0), P = O p (0, 0) and double brackets denote the subtraction of the corre- 
sponding vacuum average at the points t,x = 0. The charge (p) and spin (a) bosonic fields 
in Eq. (JT3*j) are related to the chiral fields ^(z^tq) introduced above by the simple linear 
equations 

<At(0 p ) = -j= (<Pr,i ± Vl,t T - <^,|) (14) 

(the upper sign corresponds to ip a and the lower sign denotes B p ). With the help of the 
canonical transformation Eq. (|55|) the chiral bosonic fields in Eq.(|14|) can be expressed in 
terms of noninteracting plasmonic modes q9 b; t, j (j = 1,2). For an infinitely long LL the 
propagators of these fields are (see e.g. Ref. [32|) 

(WR/L,j{ti X )¥R/L,k)) = -^ ln — • ( 15 ) 

where the velocities Sj are defined in Eqs. (|58|) and ()59|). Finally, the expression for the 
Josephson current through a side-coupled LL wire (FigHJ) takes the form 



J^\V ,Xa,A Z ^) = 

eVF Df 2 { 2(71_1) ^B(l/2, 7l - l/2)F(l/2, 7ls ; 7l ; 1 - ( S2 / Sl ) 2 ) sin (<p - g- 
2n z I V L ) s\Vp V Ail 



+ (|) 2(72_1) ^i?(l/2,7 2 -l/2)F(l/2, 72c ; 72 ;l-(^i) 2 )sin^+ , (16) 

where B(x, y) = T(x)T(y) /T(x + y) is the beta function, F(a, j3; 7; z) is the hypergeometric 
function (see e.g. Ref. |24j|). vp = (v 1F + v 2 f)/2 , 7 j = 7 j a + 7 j c (i — 1, 2) and 

w 2F sin 2 ^ cos 2 ^ 
7 is = , Tie = , 7 2, = 7i«(l ^ 2), 72c = 7 i c (l -> 2). (17) 

VlF 92 91 

Here = Sj/vjp are the correlation parameters of a two-channel LL (see Appendix I) and 
angle parameter ip is defined by Eq. (|57j) . 

By using the properties of the hypergeometric function it is easy to show that for a given 
strength of the electron-electron interaction the Josephson current jW satisfies the equations 

J^(-\ a ,A z ; V ) = J^{X a ,-A z;( p) = - J«(A«,A z ;-p) (18) 
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which describe the symmetries of electric current with respect to space and time reflections. 
In particular one can infer from Eq. (fTHj) that when both the chiral symmetry breaking 
(A a 7^ 0) and the Zeeman (Az 7^ 0) interaction are present the supercurrent can persist even 
at (p — 0. This anomalous supercurrent exists already for noninteracting electrons (Vo = 0) 
and at first we analyze Eq. in the limit of weak e-e interaction. 

For noninteracting electrons (Vo = 0,gi = <ji — 1) Eq. (jl6j) is much simplified to 

#m-*5£-(*-£) + W-(* + £)} < 19 > 

where J(°> = (Dev F /4L)(C /tt) is the critical Josephson current. We see that in the absence 
of magnetic interaction (Az = 0) the Rashba interaction in the considered geometry of SNS 
junction does not affect Josephson current at all (see also Ref. J3] where an analogous result 
was derived for a short 2D SNS junction in the presence of Rashba spin-orbit interactions). 
The interplay of the Zeeman interaction and the dispersion asymmetry in quantum wires 
results in the appearance of an anomalous (at ip = 0) Josephson current jW = J^\ip = 0) 
which it is convenient to express in terms of the asymmetry parameter A a and the magnetic 
phase xb = A Z /A L (A L = hv F /L) as 

4 l) (Aa, Xb) = ^ {(1 - Aa) sin (^-) - (1 + A a ) sin (^-) } . (20) 

As is evident from the above equation, the anomalous supercurrent J a appears only when 
both the dispersion asymmetry and the Zeeman interaction are present J a (A a = 0, Az) = 
J a (A a , A z = 0) = 0. In the limit of weak dispersion asymmetry A a 1 (a realistic case- for 
quantum wires formed in 2DEG) the Josephson current as a function of Zeeman splitting 
demonstrates a simple harmonic behavior with a slow periodically varying amplitude (beats) 

J« ~ Ji°) sh 



cos ( ^ ) . (21) 



V + Xa [A- L - tan A- L J\~~~\A L 
Now we analyze Eq. (|16J1 in the limit when the Rashba interaction is negligibly small 
(A a = 0). In this case the Josepson current through the LL wire takes the form 

JW = J^cosxbSwv , Ji 9) = R(g c ) J { "\ (22) 



where the interaction-induced renormalization factor R(g c ) (here g c 1 = J 1 + 2Vq/hTivf 
is the LL correlation parameter in the charge sector) is equivalent to the one evaluated in 
Ref. |ll I 

R(9c) ~ 7^r(l/2 + l/2, c )) F {2 ' 2' 2gZ + 2' 1 " 9c ) [l] (23) 
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In the limit of strong interaction Vq/Tlvf 3> 1 the renormalization factor is small 




) 



3/2 



( 



a 



) 




R(g c < 1) 



2 ^ V 



L 



< 1 



(24) 



and the Josephson current through a S/I/LL/I/S junction is strongly suppressed^. 

When both the electron-electron interaction and the dispersion asymmetry are strong, 
only one of the two terms in Eq. ([TH|) dominates. The corresponding critical current (for 
definiteness we assume that v\p ~ Vf/2 ^> V2f) 



is much smaller than the critical current J c in the absence of dispersion asymmetry (v%f = 
V2f)- It means that chiral symmetry breaking in quantum wires enhances the interaction- 
induced suppression of the Josephson current. 

III. DISPERSION ASYMMETRY AND QUASI-RANDOM ENERGY SPEC- 
TRUM OF PLASMONS 

In this section we evaluate the spectrum of topological excitations and plasmonic modes 
in a LL wire of the length L end-coupled to bulk superconductors (see Fig|3j). The electron 
normal backscattering at the N/S interfaces is assumed to be strong. The Josephson current 
can be calculated to the first order on junction transparency using Eq.© for the (^-dependent 
part of the ground state energy. For a finite length LL the zero modes in Eqs. Q and (fTUJl 
contribute to the energy and after some algebra we get for SE^ 2 \(p) an expression analogous 
to Eq. (JT2J) where now IIi(2)(t) is replaced by the product Hi(2)(t)Qi(2)(t). The zero mode 
contributions Qu2){t) are (J = 1,2) 



To zeroth order of perturbation theory in the barrier transparencies the electrons are confined 
to the normal region. Therefore the correlation functions in Eq. (fT2*|) have to be calculated 
with the appropriated boundary conditions. The natural boundary condition for our problem 
is the requirement that the particle current through the interfaces at x = 0, L is zero 




'IF 



(25) 



Qjir) = exp J -2tt( A (fl T - ft J + ^ (ft T + ft;) ) i exp 





(27) 
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Here the wave function ty a for the nonsymmetric electron dispersion is represented as 

# m) ~ e lk ^ x ^ Rm {x) + e- ik W x -$ mx {x). (28) 

Notice that Eqs. ()27|) and ()28|) determine more general bound ary conditions than ^ a (x = 
0, L) = usually assumed in the literature (see e.g. Ref. 2Cj). The last b.c. is the 
particular case of so called "hard wall" b.c.'s ^^\xh) = j = 1, ...,2N for a multichannel 
(N) spin-1/2 LL. They do not mix the channels and allows one to reduce the multichannel 
problem to calculations for a single channel situation with an additional summation of 
channel dependent quantities over channel quantum numbers^. In our case scattering at 
the boundaries changes the channel "index" (1 <-> 2) and the correct b.c. for "slow" fields 
*$?r(l) has to take this fact into account. The decomposition Eq. (|2"%)) holds at distances 
much larger then A^. In a general case, the wave function at the boundary is of a more 
complicated and unknown form and one may not put ^ = in order to find the relations 
between the two terms in Eq. (J28|) . In contrast, the requirement that the particle current 
through the boundary is zero is robust and its consequences hold at any distance from the 
boundary due to current conservation. 

For the bare electron spectrum without dispersion asymmetry {kip = k 2 F — kp) the 
formulated requirement is equivalent to the following boundary conditions for the chiral 
(R,L) fermionic fields (see also Ref. 

*L(aO*J*(z)l*=o,L = *L(^)^U=o,l- (29) 

The boundary conditions Eq. (J2T?|) correspond to a LL with open ends^ and result in zero 
eigenvalues of the momentum-like zero-mode operator 11^ and in quantization of harmonic 



modes (plasmons) on a ring with circumference 2L (see Ref. 22]). In this case the spectrum 



of plasmons is equidistant and the propagators take the form (j, k — 1,2) 

a \ u Sjk, 1 - exp[27r(±a: - s k t + ia)/L] 
({(pR(L)j{t,x)(p R{Ij)k )) = -^ln . (30) 

Here a is the cutoff length and Si^) are the velocities of charge and spin excitations (for 
noninter acting fermions Si = So = vp)- 

n 

Now we generalize the quantization procedure elaborated in Ref. [22J to an electron 
spectrum with dispersion asymmetry. We will assume that electron normal backscattering 
at the boundaries is not accompany by spin-flip processes. Therefore each backscattering 
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for our spectrum (FigJIJ leads to the change of the channel index ("1" <-> "2") and the 
corresponding Fermi velocity. 

It is worthwhile at first to consider the general case of boundary scattering in a two- 
channel system of noninteracting electrons confined to the interval [0,L]. The electron 
backscattering at the boundaries is described by 2 x 2 unitary symmetric matrix which 
is convenient to parameterize as follows 

„ / t i\t\\ 



S 



i\t\ 



(31) 



where r 



re 



is the intrachannel backscattering amplitude (1 <-> 1,2 <-> 2) and t is the 



interchannel backscattering (1 <-> 2) amplitude |r| 2 + |t| 2 = 1. By matching the electron 
wave functions at the boundaries x = and x = L with the help of the S-matrix Eq. (|31|) 
one easily finds the spectrum equation 



cos 



eL 
2 \v 1F 



— (— + —) 

v 2F J 



+ 5 



I 1 2 2 

r cos 



eL f 1 1 \ ' 

2 V^iF ^2F/ 



(32) 



For purely intrachannel reflection, t = 0, we get from Eq. (}3*2~j) two independent sets (j = 1, 2) 
of equidistant levels with spacing Ae^ = nhvjp/ L. In the opposite case of purely interchannel 
backscattering (r = 0) the spectrum is also equidistant 



2nTi v\fV2f ( 15 
\n H 

L Vi F + V 2 F\ 2 7T , 



n 



1,2, 



(33) 



In a general case the spectral equation (|32|) yields a set of quasi-random energy levels. 

The bozonization technique is compatible only with the two considered limiting cases: 
|r| = 1 (this was demonstrated in Ref. 221]) . and r = as we will show now. Let us start 
at first with the case on noninteracting fermions. The boundary condition Eq. (|2*7jl for 
Vif 7^ V 2F results in the equations 



v 1F ^ m (x)^ R1 (x) 



x=0,L 



v 2F ^ L2 (x)^ L2 (x) 



\x=0,L, 



V{ 2 {x)ty R1 {x)e 



i(k 1F +k 2 F)x 



x=0,L 



0. 



These equations are satisfied if 

a 2 L 2 



VlF 
ViF 



1 1 1 

L\ L 2 L 



F 27T V 1F V 2F 

sZ = — n 



L Vif + V 2 F 



n 



1,2,... 



(34) 
(35) 

(36) 
(37) 
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and 



[$ L *(x,t) + $ R(7 (x,t)]\ x=0jL = ^n a , a=U, (38) 



where and are integers. Eq. (|3*5j) in its turn is satisfied for topological sector with 
quantum numbers (<^Ro- + ( pLa)/\/^ — n a, iX = and the harmonic modes ^R(L)u{x,t) 
which obey the relations 

(fRa{x,t)\ x=0)L = -<PLa(x,t)\ x =0,L- (39) 

^From Eqs. (fTT|) . (j3Tlj) and (J3~9"j) one easily gets the plasmon spectrum 

e n = ■ n (40) 

L si + s 2 

(si t 2 are the plasmon velocities, which coincide with the Fermi velocities for noninteracting 
fermions) and the desired correlation functions (j, k = 1,2) 

// , ,x u S jk 1 - exp[z27r(±x - s k t + ia k )/L k ] 
(( W)j (x lt ) W)k )) = -^ln 2nak/Lk 1, (41) 

where the effective quantization lengths Lj according to Eq. (}3l)j) are 

T Vi F + V 2 F T /<)rA 

L H 2) = L (42) 

V2(1)F 

In the limit v\f = V2F Eqs. (|47H) - (|4*2"|) reproduce the plasmon spectrum and the correlation 
functions of a single channel LL with open ends^. 

Now we are ready to consider the effects of interaction. For a single-mode LL (or for 
a multichannel LL, provided the backscattering is allowed only to its own channel) the 
boundary condition Eq. (|3*9*j) for harmonic modes holds also for interacting fermions as one 
can check using a Bogoliubov-like transformation which diagonalizes the LL Hamiltonian. 
Hence in the presence of interaction one can still use the same correlation functions as 
for noninteracting fermions with the only difference that the velocities are renormalized by 
interaction. 

This is not the case for our problem. With the help of exact transformations (see Eq. (jSU)) 
which diagonalize the 2-channel LL Hamiltonian^ one can show that if the chiral bosonic 
fields satisfy Eq. (|3U|) . the diagonalized ones <fR(L)j are connected at the boundaries by the 
effective "scattering matrix" S e 

(p Rj {x = ^L) = Y.S%(p Lk {x = ^L) , S» = lf~ ^, (43) 
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where the coefficients A, B are 



A — — [cos 2ip - sinh($i - i? 3 ) sin2'0]" 1 [sinh(^i - i? 2 ) cosh($i - # 3 ) (44) 
+ cos 2^ cosh(#i — i? 2 ) sinh(i?i — i? 3 ) + sin 2^/> cosh(#i — i? 2 )], 

B = - [cos2^-sinh(7?i - ^ 3 ) sin 2^] - 1 [cosh (^! - t? 2 ) cosh(i?i - i? 3 ) (45) 
+ cos 2ip sinh(i?i — $ 2 ) sinh(i?i — $ 3 ) + sin 2ip sinh($i — i? 2 )] 

and the "rotation angles" $1 (Z = 1, ...4) and ^ are defined in Appendix I (see Eqs. ()56|) and 
f)57p ). One can check after some algebra that the coefficients A and B satisfy the simple 
relation B 2 — A 2 = 1, which makes the S-matrix in Eq. (|43|) unitary. This observation allows 
us to use the scattering matrix formalism when evaluating the energy spectrum of plasmons. 
Notice that in the parametrization Eq. we have r = iA/B, \t\ = 1/B, 5 = tt/2. 

For monochromatic bosonic fields with amplitudes &_r(l)j the scattering at the boundaries 
x = and x = L are determined by the equations 

^ = 0: b Rj = J2S] k b Lk , x = L: e~^b Li = £ S° k e ia *b Rk , (46) 

k=l k=l 

where the phases ctj = sL/sj and Sj are the plasmon velocities (see Eqs. (JH%J) . (JHnj) ). ^From 
the above set of linear equations one easily finds the spectrum equation for plasmons 

















[t 






= Rsm 2 


[f 





















(47) 



where R = (A/B) 2 < 1 is the effective backscattering coefficient for plasmons. It depends 
both on the dimensionless interaction strength k = Vo/kTi{vif + u 2 f) and on the dispersion 
asymmtry parameter A a . Notice that the spectral equation (|47jl is the special case of Eq. (132}) 
for 5 = 5 r = 7r/2. 

The derived spectral equation has simple exact analytic solutions in two limiting cases: 
(i) noninteracting fermions and, (ii) when dispersion asymmetry is absent, V\p = V2F — Vf- 
For noninteracting particles (Vo = in Eq. ( flj) ) the "rotation angles" are ip = 0, $1 = $4 = 
(see Appendix I) and the velocities S\ = Vip, s 2 = v 2F- Then A = 0, B = — 1 and the 
effective backscattering coefficient R = 0. Eq. (}4Tj) in this limit reproduces the equidistant 
spectrum of "noninteracting" plasmons, Eq. ()33j) . For interacting fermions in the absence 
of dispersion asymmetry the "rotation angles" are $1 = t? 3 ,i? 2 = , &4,cos2ip = 0. In this 
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limit Si = s,s 2 = v F and R — 1. The corresponding plasmon energies e£> = nsn/L, = 
nvFn/L, (n = 1,2, ...) represent the standard excitations in the charge and spin sector of a 
finite length (L) Luttinger liquid with open ends^. 

For a general case Eq. ()47|) has to be solved numerically and the plasmon spectrum 
represents a set of quasi-random energy levels. The plasmonic energies can not be separated 
into two independent sets of levels - one for charge density excitations, another for spin 
density excitations. It means that the considered boundary conditions strongly mix the 
charge and spin excitations and the phenomena of charge-spin separation, well known in a 
LL, generally speaking, can disappear when both spin-orbit interaction and finite size effects 
are present. 



IV. JOSEPHSON CURRENT THROUGH A FINITE-LENGTH LL WIRE WITH 
DISPERSION ASYMMETRY 

It is clear from a physical point of view that the effects of a dispersion asymmetry in 
the bare electron spectrum have to be most significant in the quantum dynamics of non- 
interacting electrons. In this case the mismatch in Fermi velocities when an electron is 
backscattered at the boundaries leads to an intricate interference pattern. The more strongly 
particles interact, the less important are the effects of dispersion asymmetry. For instance, 
in the limiting case of a ID Wigner crystal (strong repulsive long-range interactions) it is 
hard to imagine any interference produced by the quantum dynamics of plasmons in two 
Wigner crystals pinned by structural imperfections at the boundaries. So in our problem it 
is reasonable to expect the restoration of the regular plasmon spectrum and the spin-charge 
separation in the limit of strong interaction. 

For strongly interacting electrons k = Vq/itTi{v\f + v 2 f) 1 and v\p ~ v 2 f (i-e. for 
week or moderate dispersion asymmetry) the coefficient R (intra-channel backscattering 
probability) in Eq. (jTTj) takes the form 

r^L-L W*-*® (48) 

R ~ #>• 2^ ' (48) 

We see that when k ^> 1 the backscattering is always an intrachannel process (R ~ 1 with 
a high accuracy) and the spin-charge separation and the e quid istant character of plasmon 
spectra are indeed restored. This observation (see also Ref. |23|) justifies for strongly inter- 
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acting multichannel (j = 1,...N) LL the boundary conditions = 0, L) — 0) usually 

postulated in the literature (see e.g. Ref. [2(|) for arbitrary interaction. 

It is straightforward to evaluate the Josephson current using the exact plasmon spectrum 
for R = 1 (i.e. when A a = 0) and the propagators Eq. (jHO)) . The result for zero Zeeman 
splitting (Az = 0) is 



where J(°> = (Dev F /4L){C fir) is the critical current through a ID SNS junction and the 
interaction- induced renormalization factor Rf(g c ) is 



Here F(a, (3; 7; z) is the hypergeometric function and g c is the LL correlation parameter (see 
Eq. (|22p). For noninteracting electrons Rf(g c — 1) = 1 and our formula has to reproduce the 
known expression for the Josepson current through a ID SNS junction (see e.g. Ref. 
From this comparison one finds C = 7T. 

As we have already mentioned in this section, R — ► 1 in the limit of strong interaction 
k ^> 1. This observation allows us to evaluate the correlation functions and the Josephson 
current for strongly interacting electrons with dispersion asymmetry. The Josephson current 
is described by Eqs. ()49|) and (joT?|) after the replacement g c — > n~ l l 2 and in the limit k ^> 1. 
The renormalization factor now takes the form 



The formulae (J49)) and (|51|) show that in the considered limit the Josephson current dos not 
depend on the parameter A a of dispersion asymmetry. By comparing Eq. (|HT|) and Eq. (|25j) 
we see that the interaction suppresses supercurrent more strongly in a long end-coupled 
quantum wire then in a side-coupled one, which is in agreement with intuition. 

Dispersion asymmetry affects the supercurrent of weakly interacting electrons. The influ- 
ence, however, numerically is not strong even for the most favourable case of noninteracting 
particles. With the help of the correlation functions (f4"T|) it is easy to calculate the Josephson 
current of noninteracting electrons with dispersion asymmetry 



jW(g c ;<p) = jWR f (g c ) simp, 



(49) 




(50) 




(51) 




(52) 
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Here Jf^ is the critical current in the absence of dispersion asymmetry (see Eq. ([49))) and 
R(X a ) is the renormalization factor induced by the asymmetry of electron dispersion 

^(i-y Ji + (,76-da; 

sm(7rA a ) A a -> 1 

We see from Eq. that the dispersion asymmetry always slightly enhances the critical 
current. The analysis of the Josephson current in a ID SQWS junction in the presence 



2l\ using the Andreev level approach. It 



of dispersion asymmetry was performed in Ref. 
was shown that the observed enhancement of the Josephson current is due to less perfect 
cancellations (different Fermi velocities) of the partial supercurrents carried by adjacent 
Andreev levels. 



V. THE RASHBA EFFECT, CHIRAL ELECTRONS IN ID QUANTUM WIRES 
AND THE JOSEPHSON CURRENT IN SLLS JUNCTION 

Now we consider the limit of strong Rashba interaction. In this case the electrons in 
a quasi-lD quantum wire behave like truly chiral particles, that is the spin polarization 
of an electron irrespective of its subband index is determined by the direction of elec- 
tron motion along the wire — right-moving and left-moving electrons have opposite spin 
projections^. We will assume for defmiteness (it depends on the sign of Rashba coupling) 
that "R"-electrons are "down" -polarized and "L"-electrons are "up" -polarized (see FigEl 
dashed lines for spin polarizations). We have already seen in section II, that the left/right 
symmetry breaking in the presence of the Zeeman interaction results in the appearance 
of an anomalous Josephson current. Physically it means that when the spin projection is 
correlated with the direction of motion (left, right), the magnetic field, via the Zeeman in- 
teraction, induces partial Josephson currents (for each subband "1" and "2") even if the 
superconducting phase difference in the SNS junction is zero. For the spin alignments as- 



sumed in Refs. |l6lll7l ] the subbands contribute to the Josephson current with opposite signs 
and therefore the anomalous supercurrent vanishes for symmetrical spectrum v if = V2F- As 
we see, the electron dispersion asymmetry is indispensible property to get anomalous Josep- 
son current for a weak Rashba interaction. In the limit of strong s-o interaction when all 
right (left) moving particles have parallel spins, the contributions of subbands have the same 
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sign and the existence of electron dispersion asymmetry ceases to be crucial in appearence 
of anomalous (at ip = 0) Josepson current. 

What is more important are the spin-flip processes which may take place in the transition 
regions between the 2D or 3D electron reservoirs (superconducting leads in our case) and 
the ID quantum wire with a pronounced Rashba effect. The electrons in the reservoirs have 
two possible spin states, while deep inside the wire, where the s-o interaction is strong, the 
electron spins have to be aligned according to the above discussed prescription. So particles 
with the "wrong" spin projection should turn their spins toward the "right" direction. 

One can imagine two different types of transition regions. In the case when s-o interaction 
is changed abruptly at the lead/wire interfaces, the spin-flips induced by the Rashba inter- 
action will be accom pan ied by normal electron backscattering. Such non-adiabatic contacts 
were studied in Ref. |l8| when evaluating normal electron transport through a ID quantum 
wire with strong Rashba interaction attached to leads where the s-o interaction was assumed 
to be negligibly small. In this model the transparency of the junction strongly depends on 
the spin-orbit coupling and normal electron transport is suppressed even for perfect contacts. 

Here we will assume that the Rashba interaction in the QW is switched on smoothly and 
that therefore the backscattering of the electron at the boundary and the rotation of its spin 
induced by the Rashba interaction the in quasi-lD quantum wire are spatially separated. 
For instance, the left-moving electron (spin-" up") at first is normally backscattered at the 
left interface keeping the spin-" up" and then after travelling some length <C I <C L its 
spin becomes "down" -polarized. In this model, the Rashba interaction does not lead to 
additional electron backscattering at the interfaces and does not suppress the supercurrent. 

It is also convenient for calculations to assume that the electron-electron interaction in 
the ID QW is strong. As was shown in section IV, the interchannel electron backscattering 
at the I/LL interfaces is suppressed in the limit of st rong repulsive interaction and one can 
use a simple quantization procedure proposed in Ref. 22] to evaluate the plasmon spectrum 
and the correlation functions. After straightforward calculations the desired expression for 
the Josephson current takes the form 



J (R \ V ) ~ J®R f sin 



■ A z f 1 1 



cos 



A z f 1 I 



2 VA L1 A L2 



(54) 



2 VA L1 A L2 , 

where the interaction- induced renormalization coefficient Rf is determined by Eq. 1)51)1 . As 
was already evident from physical considerations, the anomalous supercurrent J^ R \cp = 0) 
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in the limit of strong Rashba interaction is induced by a magnetic field (Az 7^ 0) even in 
the absence of electron dispersion asymmetry. We see from Eq. (I54j) that the dependence 
of the supercurrent on magnetic field is absolutely different for chiral and normal electrons. 
In particular the critical current for symmetric electron spectrum (vip — v 2f) in the case 
of chiral electrons does not at all depend on the Zeeman splitting, while in the ordinary 
situation one gets a periodic dependence. 

VI. CONCLUSION 

The problem we have studied allows one to consider the interplay of proximity-induced 
superconductivity, the Rashba, Zeeman and Coulomb interactions on the transport proper- 
ties of quasi- ID quantum wires. We have shown that the Rashba and Zeeman effects strongly 
influence the supercurrent. In particular, the Rashba effect in quantum wires results in a 
strong correlation between electron spin polarization and the direction of electron motion. 
In other words a strong Rashba interaction creates chiral particles in the ID electron system. 
The influence of a magnetic field via the Zeeman interaction on chiral particles leads to the 
appearance of a net electric current in the wire. When the leads that the quantum wire is 
attached to are superconducting this current persists even at zero phase difference across the 
junction. The effect exists already for noninteracting particles and it is strongly sensitive to 
the electron dispersion asymmetry for weak Rashba coupling and ceases to depend on the 
asymmetry parameter in the regime of strong Rashba interactions. 

It is well know n 10 i 26 that the Josephson current in a perfectly transmitting junction (i.e. 
without normal electron backscattering) is not influenced by the Coulomb interaction. In 
contrast, any potential barrier inside the normal region which induces electron backscat- 
tering is renormalized (upwards) by the repulsive interaction (the Kane-Fisher effect) 
and therefore strongly suppresses the supercurrent through a (poorly transmitting) SIL- 
LIS junctio n 10 ! 11 ! 20 ' 26 . We have shown that the electron dispersion asymmetry, which is 
induced by the Rashba interaction in quasi-lD quntum wire a 16 ! 17 is significant for the super- 
conducting properties of a LL junction only for weak or moderate Coulomb interaction. In 
this case the interplay of interaction and dispersion asymmetry leads to an intricate interfer- 
ence pattern in the plasmon quantum dynamics of a finite length two-channel LL and makes 
the plasmon spectrum quasi-random. Strong Coulomb interactions suppress this kind of 
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quantum behavior and restores a regular (equidistant) plasmon spectrum. The tendency of 
strong Coulomb interactions to suppress quantum interference can be traced in different ID 
electronic systems, for instance, in a LL double barrier (absence of resonant tunnelling for 
a strong repulsive interaction^) or in mesoscopic coupled rings (ordering effect of Coulomb 
interaction on persistent current oscillations^). 
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Appendix I 

The canonical pseudo-orthogonal transformations, which diagonalize the Luttinger liquid 
Hamiltonian 



are^ 3 - 



f Pr\\ 

Phi 
PRI 
Vplt / 



/ cosh $1 cos if) sinh $1 cos if> — cosh {} 2 sin if> — sinh d 2 sin ip \ 

sinh cos ^ cosh $1 cos if) — sinh $2 sin if) — cosh $2 sin if) 

cosh $3 sin ^ sinh $3 sin if) cosh $4 cos if) sinh $4 cos if) 

V sinh $3 sin if) cosh $3 sin if) sinh $4 cos r/> cosh $ 4 cos ^ / 



( Pri\ 
Plx 

PR2 
\pL2J 



(55) 



where the "rotation angles" dj and if) are expressed in terms of the Fermi velocities V\ F , v 2F 
and the interaction strength V by the following equations 



01 



1 



■ln^i, $ 2 



In 



f V\F 
\V 2F 



92 



a 



In 



f V2F 
\VlF 



9i 



0, 



1 



In g 2 , 



tan 2-0 



2V y/Vl F V 2 F 



(56) 



(57) 



(f IF - V 2F ) [Vq + 7T^(f IF + «2f)] 

Here gr 3 - = Vj F /sj (j = 1, 2) are the correlation parameters of a 2-channel LL and the plasmon 
velocities are 

1/2 



Si = Vip { cos V + 



\v 1F J 



sin -0 + 



Vn 



irhvxp 
Si(lf> <-> -ip,v 1F 



cosip + \ sinip 

V % 



(58) 



s 2 = <-> <-> U2f ,). (59) 

For noninteracting electrons, Vq = 0, the correlation parameters are #i = #2 = 1 and, 
according to Eqs. (jST?)) and ({37}) $1 = $4 = 0, if> = 0. In the limit V\ F = v 2F = v F Eqs. (f5T)))- 
(jSTJj) reproduce the well-known expressions for the correlation parameters of a spin-1/2 LL 



01=^3 



02 = 04 = 0, £ C 



2Vp \ 

7rhv F J 



-1/2 



(60) 
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